Disordered Totally Asymmetric Simple Exclusion Process: Exact Results 
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We study the effect of quenched spatial disorder on the current-carrying steady states of the totally 
asymmetric simple exclusion process with spatially disordered jump rates. The exact analytical 
expressions for the steady-state weights, and the current are found for this model in one dimension. 
We demonstrate how these solutions can be exploited to study analytically the exact symmetries of 
the system. In particular, we prove that the magnitude of the steady-state current is left invariant 
when the direction of all the allowed particle jumps are reversed. Or equivalently, we prove that for 
any given filling and disorder configuration, particle-hole transformation is an exact symmetry that 
leaves the steady-state current invariant. This non-tri vial symmetry was recently demonstrated in 
numerical simulations by Tripathy & Burma (preprint cond-mat/9711302). 
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Driven diffusive systems have been a subject of exten- 
sive studies in recent years [0 . A key ingredient of these 
systems is the presence of a driving field due to which 
the underlying dynamics do not generically obey detailed 
balance, leading to steady states with non- vanishing cur- 
rents. In the absence of a unifying theme that encom- 
passes our understanding of nonequilibrium phenomena, 
these systems offer a relatively simple framework within 
which systems far from thermal equilibrium may be stud- 
ied. 

Driven diffusive systems further provide the possibil- 
ity of studying the intriguing interplay of disorder, inter- 
actions and drive for a wide range of parameter values 
like the degree of disorder, filling and drive strength [Q. 
These expectations are not unfounded: (i) Driven diffu- 
sive systems in the absence of disorder have been stud- 
ied extensively and have revealed basic differences be- 
tween equilibrium and nonequilibrium systems fl],[|,|j. 
For instance, it has been demonstrated that sponta- 
neous symmetry breaking (SSB) |||(| and phase sepa- 
ration can occur in one dimensional asymmetric exclu- 
sion processes Q. In contrast it is well known that one 
dimensional systems in thermal equilibrium with short 
range interactions do not exhibit phenomena like SSB 
and phase separation S; (ii) Also, systems with disorder 
and drive but no interactions between particles are well 
studied and understood 0. 

A number of physical situations involving flow in ran- 
dom media require an understanding of disorder driven 
diffusive systems of interacting particles [|l0) . Most of the 
understanding gained in these system have been largely 
based on numerical simulations. Analytical character- 
ization in terms of the exact steady state measures in 
systems without translational invariance have been found 
only in the case of the disordered drop-push model [jllj . 

An interesting class of models for which there exists no 



analytical characterization of the steady-state weights or 
the steady-state current is the disordered totally asym- 
metric simple exclusion process (DTASEP). Even the sin- 
gle disorder case has not been amenable to an analytic 
treatment. Extensive numerical and mean-field studies 
have been done and a number of interesting features have 
been highlighted JH|-[lf]. 

In this work we derive, for the first time, formally exact 
solutions for the steady-state weights, and current for the 
DTASEP model in one dimension. These solutions are 
valid for arbitrary disorder and particle fillings. We fur- 
ther show how these solutions can be exploited to study 
the exact symmetries of the model by studying the trans- 
formations that leave the magnitude of the steady-state 
current invariant. In particular we show that the model 
possesses an exact particle-hole symmetry for any given 
disorder realization; a result that was observed in numer- 
ical simulations by Tripathy & Barma p| . In addition 
we prove the "obvious" result that the steady-state cur- 
rent is a constant across each bond. Although this result 
is intuitively obvious, given that the dynamics conserves 
the number of particles, to the best of the authors knowl- 
edge no explicit proof of this statement for a disordered 
lattice model exists. 

Definition of the model: The DTASEP model is de- 
fined on a ID lattice of length L with periodic boundary 
conditions. Each site can hold either 1 or particle. 
Each bond + of the lattice is assigned a quenched 
random rate on chosen independently from some chosen 
probability distribution. The evolution is governed by 
random sequential dynamics defined as follows: in a time 
interval dt the particle attempts to hop, with probability 
ctidt to its neighboring site i + 1. Wc consider the case 
in which the jumps are allowed only in one direction (to 
the right) and is the same for all bonds. In addition, the 
move is completed if and only if site % + 1 is unoccupied. 
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The time averaged steady-state current Ji in the bond 
(i, i + 1) is given by: 



Ji = ai{rii(l 



(1) 



where rij counts the number of particles at site i. 

Since for every particle hopping to the right, a "hole" 
jumps to the left, implies interchanging all the particles 
for holes (charge conjugation - C) and reversing the di- 
rection of hopping (time reversal - T) leaves the steady 
state current, up to a sign, the same. Hence the steady- 
state current is symmetric under a CT transformation. 
This symmetry is valid in general: in any dimension; in 
the presence of disorder; and in the case when the parti- 
cle can hop along any direction with finite probabilities 
(DASEP). 

In Eq. ([l]), if all the a^s are put equal to the same con- 
stant a, then it follows that the steady state current Jo, 
which is assumed to be the same on each bond and hence 
Jo = L zZi Jii i s invariant under m — > 1 — n% for all 
i. Hence for the clean totally asymmetric simple exclu- 
sion process in ID (TASEP), charge conjugation (C) in 
itself leaves the steady-state current invariant. Since CT 
symmetry is always true, it follows that for the TASEP 
model in ID both charge conjugation and time reversal 
(T) symmetries individually holds. 

The surprising observation made in the numerical sim- 
ulations in Ref. | fl4| was that, charge conjugation symme- 
try was found to hold even in the presence of disorder in 
the DTASEP model in ID. This symmetry as is obvious 
from Eq. ([!]) is not expected to hold when all the ai 's are 
taken to be random. Convincing numerical evidence has 
been offered, although a general proof of the validity of 
this symmetry for all fillings and disorder configurations 
has not been obtained. 

We derive below for this model a formally exact ex- 
pression for the steady-state probability density and the 
current Jo. Using these expressions we show that the 
observed reflection symmetry holds in general for all fill- 
ings and disorder configuration which in turn implies the 
symmetry of the steady-state current under charge con- 
jugation. 

For L sites with N particles, the number of configu- 
rations M = L C N . Here L C N = L\/(L - N)\N\ is the 
number of ways N particles can be distributed amongst 
L sites with a maximum of only one particle per site. The 
dynamics of these configurations, for a given realization 
of the quenched bond variables R = {ai, 0:2, • • • , <Xl}j 
are given by the rate equation: 



dP m 

dt 



T(n — > m)P n — T(m — > n')P n 



(2) 



where T(m — » n) gives the transition probability from 
configuration m — > n. This expression can be conve- 
niently expressed in matrix notation as: 



dP(t) 
dt 



= WP{t) 



(3) 



with W mn = T(n — ► m) and W mm = - J2„> T(m —> n'). 
Since 1 • P(t) = ^2 m Pm{t) = 1 V t (normalization), it 
follows that 1 • W = 0. This implies that W has a zero 
eigen- value with left eigen- vector 1 . Hence there also ex- 
ists a right eigen- vector P such that WV = 0. This de- 
fines the steady-state solution of the problem for a given 
choice of the disorder R. 

The steady-state of the problem considered above is 
characterized by a uniform current Jo across each bond. 
Given the solution P, the current across, say, site i and 
i + 1, is given as Y]Pm = Jo, where the sum is only 
over the set of configurations with a particle on site i and 
a hole on site i + 1 (see Eq. 0) . 

The above sum for all sites i can be combined and con- 
veniently written in matrix notation by defining a ma- 
trix G such that G ■ P = Jg(l/ai, l/a2, • • • , l/ai J ) T = 
Jo(l/a) T . The entries of G for a given row i (corre- 
sponding to site i) has I corresponding to the configu- 
rations with a particle at site i and a hole at site i + 1 
and otherwise. The order of the G matrix is therefore 
[L x M) since there are L sites and M configurations. 
The G matrix can be expanded to an [M x M) matrix 
by augmenting an (M — L) x M zero matrix to G, such 
that: 



G L 

Or, 



Jo 



1/q 




JoV (4) 



where for convenience we have defined new variables 
Q = (G,0) T and V = (l/a,0) T . Hence the three equa- 
tions that determine the steady-state distribution and 
current are: 



W-P = 
Q ■ P = J V 
1 • P = 1 



(5a) 
(5b) 
(5c) 



If we assume that the row-rank(M / ) = M — 1 and 
not lesser, then there exists a one-parameter solution to 
Eq. (|5a|). The most general form of the solution will be 
of the form /3P, where /3 is the arbitrary parameter. The 
solution when substituted into Eq. ([5^) allows the f3 fac- 
tor to be absorbed in the definition of Jo by re-scaling 
Jo - ¥ Jo/0- Hence, the most general solution P = P(Jo) 
will have Jo appearing as the only free parameter. The 
value of Jo is fixed by the normalization condition given 
in Eq. ([k]). Hence a unique P and Jo solves Eq. (0). 

It should be noted that both Eq. |5a| and Eq. |3t| can- 
not be inverted to obtain a solution for P and Jo as 
neither W nor Q are invertible matrices. Although, 
W and Q are independently non-invertible, their sum 
Q + W is invertible. The proof of this is simple. Suppose 
row-rank(Q + W) = r < M. Then the equation: 
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(Q + W) ■ P = J V 



(6) 



obtained by adding Eqs. (|5a|) and ([5t]) would have a gen- 
eral solution P( Jo, Pi, 02, ■ ■ ■ , /?M-r), i-e., in addition to 
the free parameter Jo , there will be M — r other free pa- 
rameters. However, it was shown earlier that the most 
general solution of P = P(Jo) can have only one free pa- 
rameter Jo, that can be fixed by normalization. It follows 
that r = M which implies that (Q + W) is invertible. 

Using Eq. (§), we obtain P(J ) = JoiQ+W^-V. The 
steady-state current Jo is obtained using the normaliza- 
tion condition in Eq. (|5c|) giving, 1/J = l-^Q+W)^ 1 -V. 
Hence the final solutions are given as: 



P 

Jo = 



(Q + WQ- 1 • V 
1 • (Q + W)- 1 ■ V 
1 



1-iQ + W)- 1 -V 
Note that in general a linear combination: 
(Q + r)W) ■ P = J V 



(7a) 
(7b) 

(8) 



could have been used. This is equivalent to the equation 
(Q + W)-P + (t)- l)W • P = J V. Since W • P = we 
get back Eq. (g) implying that P will be independent of 
•q as it should be. 

In arriving at the above solutions, we had assumed 
that in steady-state, the current across each bond was 
the same constant Jo- The arguments presented above 
clearly indicate that there would be no unique solution 
had we started by assuming that the currents were dif- 
ferent. For suppose the currents across the bonds were 
{ Ji, J2, • • ■ > Jr,}) then the most general solution would 
be of the form P( Ji, J2, • • • , Jl). Since there are L un- 
knowns but only one normalization condition, a unique 
solution would not exist. We therefore conclude that the 
currents in steady state are a constant across each bond. 

We now use these solutions to prove the reflection sym- 
metry of the current observed in Ref. JL4]. It was, as 
stated earlier, observed in Monte Carlo simulations that, 
for a given realization of the quenched bond variables 
R = {ai, a2, • • • , the steady state-current was the 
same whether the particles were allowed to jump to the 
right or to the left. 

We denote the set of quenched random bond variables 
when the particles are allowed to jump to the left as 
R = {ai, «2, • • • , Q^l}- In R the magnitude and the or- 
dering of the ai's are the same as in R although the direc- 
tion of hopping has been reversed. Let P be the steady- 
state weights and W and Q the corresponding matrices 
such that W ■ P = and Q ■ P =_ J V. The V vec- 
tor in both cases are the same and Jo is the appropriate 
current. 

The question we ask here is: Does there exist an in- 
vertible matrix S, such that: 



Q-S = Q 
S 1 P = P 



(9a) 
(9b) 



The existence of such an S would imply that Jo = Jo, 
since: 



'■P = Q-SS- 



QP 



(10) 



We first derive the properties that S should satisfy. 
Since each row of G contains as many l's as the number 
of configuration for which there exists a particle at site 
i and a hole at site i + 1, we have G ■ 1 = L ~ 2 CV-il- 
The same is true of G, although the exact location of 
each of the l's will differ since the direction of hopping 
is reversed. Going back to Eq. (pa]), we get: 



Q-S 



— Oat- 



(11) 



All the conditions in the above equation can be satisfied 
if we take S ■ 1 = 1 

We now analyze Eq. (^) to put further constraints 
on S. Since both P and P are normalized, we obtain 
1P = 1- 5 P = 1 This equation can be satisfied if we 
take 1 • S = 1 + A, with the additional constraint that 
A • P = 0. We observe that the equation is satisfied only 
up to a vector A which is 'orthogonal' to P. 

Acting by W on both sides of Eq. (pb|), gives: 



W P = WS • P = 



(12) 



This equation can be satisfied if we take WS — ±W. The 
sign can be fixed by observing that WS ■ 1 — W ■ 1 — 
±W ■ 1. Where we have used the property that S ■ 1 = 1. 
Since the direction of hopping has been reversed in R, 
the outgoing configurations w.r.t. a given initial particle 
configuration are mapped to the incoming configurations 
in R and vice versa, we obtain W ■ 1 = —W ■ 1 (refer to 
Eqs. d&l)). 

Thus the two equations that S should satisfy are: 



QS = Q 
WS = -W 



(13a) 
(13b) 



An important feature to be noted is that it is not a simi- 
larity transformation that relates the two systems defined 
by W and W. By adding the two equations, we obtain: 



(Q + W)- 1 -(Q-W) 



(14) 



Choosing r\ = -1 in Eq. (§) and P = J Q (Q - IFT 1 • V 
with 1/J = 1 • (Q - Wy 1 - V from Eq. (|), it can be 
shown that S indeed satisfies Eq. (||) . This completes our 
proof that the steady-state currents for R and R are the 
same. 

In conclusion, following an algebraic procedure, we 
have derived formally exact expressions for the steady- 
state probability density P and the current Jq for the 
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totally asymmetric disordered simple exclusion process 
in one dimension. We explicitly show how these solu- 
tions can be exploited to study exact symmetries that 
leave the steady-state current invariant. In particular we 
have demonstrated that the steady-state current, up to 
a sign, is the same when the direction of allowed jumps 
is reversed. We also explicitly show that in steady-state 
the current across each bond is the same constant, a re- 
sult intuitively obvious but lacking proof for disordered 
lattice models. Although not showed explicitly here, the 
choice of S in Eq. (||) is a special case of the more gen- 
eral requirement that Q ■ P = Q ■ P, which guarantees 
the currents to be the same under a symmetry opera- 
tion. Hence, we believe that the methodology we have 
developed has wider applicability in studying and un- 
derstanding the symmetries of more general models not 
amenable to direct analytical studies. 

On completion of this work, we were made aware of 
an alternate proof of the time reversal symmetry of the 



DTASEP model given in Ref. |15| which in contrast to 
our algebraic proof is probabilistic in character. 
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